High-energy break-up of 6 Li as a tool to study the Big-Bang nucleosynthesis reaction 

2 H(a, 7 ) 6 Li 
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The recently claimed observations of non- negligible amounts of 6 Li in old halo stars have renewed 
interest in the Big-Bang Nucleosynthesis (BBN) of 6 Li. One important ingredient in the predicted 
BBN abundance of 6 Li is the low-energy 2 H(a,7) 6 Li cross section. Up to now, the only available 
experimental result for this cross section showed an almost constant astrophysical S'-factor below 
400 keV, contrary to theoretical expectations. We report on a new measurement of the 2 H(a,7) 6 Li 
reaction using the break-up of 6 Li at 150 A MeV. Even though we cannot separate experimentally the 
Coulomb contribution from the nuclear one, we find clear evidence for Coulomb-nuclear interference 
by analyzing the scattering-angular distributions. This is in-line with our theoretical description 
which indicates a drop of the S'24-factor at low energies as predicted also by most other models. 
Consequently, we find even lower upper limits for the calculated primordial 6 Li abundance than 
before. 



PACS numbers: 25.40.Lw,25.60.-t, 25.70.De, 26.65.+t 



I. INTRODUCTION 

The Big-Bang model of the Universe is mainly sup- 
ported by three observational evidences: the expan- 
sion of the Universe, the Cosmic Microwave Background 
(CMB), and the primordial or Big-Bang Nucleosynthe- 
sis (BBN) of light nuclei like 2 H, 3 - 4 He, and 7 Li with 
one free parameter, the baryonic density of the Universe, 
fib- Recently, a precise value for this free parameter 
(O b h 2 =0.02260±0.00053) has been deduced from the a- 
nalysis of the anisotropies in the CMB as observed by 
the WMAP satellite fl. Using the best available nuclear 
reaction rates, this now allows to make precise predic- 
tions for the primordial abundances of these light nuclei. 
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A comparison between calculations and observations re- 
veals good agreement for helium, excellent agreement for 
deuterium, and a discrepancy (by a factor of w 4) for 
7 Li 0-0] • Possible reasons for this discrepancy for 7 Li 
have recently been discussed e.g. by Spite and Spite @. 

In 2006, high-resolution observations of Li absorption 
lines in some very old halo stars have led the authors 
to claim evidence for large primordial abundances also 
of the weakly-bound isotope 6 Li Q. The 6 Li/ 7 Li ra- 
tios of ~ 5 x 10~ 2 were found to be about three orders 
of magnitude larger than the BBN-calculated value of 
6 Li/ 7 Li ~ 10~ 5 . This observation has triggered many 
studies to resolve the discrepancy either by considering 
an early 6 Li formation in primitive dwarf galaxies at high 
redshift in a hierarchical-structure formation context [7| , 
in situ by solar- like flares [8[ , or in terms of physics be- 
yond the standard model of particle physics (see, e.g. , 
Refs. @-[Hj])- More recently, however, Cayrel et ah [HJ 
and Steffen et ah [HI have pointed out that line asym- 
metries similar to those created by a 6 Li blend could also 
be produced by convective Doppler shifts in stellar at- 
mospheres. Similarly, a recent study of Garcia Perez et 
ah [14| could not claim any significant detection of 6 Li in 
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metal-poor stars. So, presently the debate is open. More 
stellar observations are required to solve this question 
(see Asplund and Lind (l5j). 

Predictions for the production of 6 Li in BBN require 
precise measurements of the 2 H(a,7) 6 Li reaction rate, the 
key production mechanism. In BBN, this reaction occurs 
at energies in the range 50 keV < E cm < 400 keV [IB] ■ At 
higher energies, this reaction has been studied carefully 
in direct kinematics: at energies above 1 MeV by Robert- 
son et al. [13], and by Mohr et al. [H| in the energy range 
around the dominant 3 + resonance at E cm = 0.711 MeV. 
At BBN energies, however, direct measurements are dif- 
ficult due to extremely low cross sections (about 29 pb 
at E cm = 100 keV). An attempt by Cecil et al. [li| at 
E cm = 0.053 MeV yielded only an upper limit for the S- 
factor of 2.0 x 10" 7 MeV b which is more than an order of 
magnitude higher than present estimates. A straightfor- 
ward solution to overcome this problem is offered by the 
indirect method of Coulomb dissociation (CD). As will 
be shown below, the dominant multipolarity involved is 
E2. When bombarding a high-Z target like 208 Pb with a 
medium-energy (« 150 A MeV) 6 Li beam, an intense flux 
of virtual E2 photons is created that dissociates 6 Li into 
2 H and a with a greatly enhanced cross section. From 
the energy-differential CD cross section, the radiative- 
capture one can be calculated easily [2Qj provided that 
the multipolarity of the respective transition is known 
and that higher-order electromagnetic or nuclear contri- 
butions can be either ignored or taken into account quan- 
titatively. 

Kiener et al. [2l[ have investigated the 2 H(a,7) 6 Li re- 
action by means of the CD method employing 26 A MeV 
6 Li projectiles breaking up into D+a in the Coulomb field 
of a 208 Pb nucleus. Referring to a theoretical paper by 
Shyam et al. [22T ]. Kiener et al. have claimed that their 
measurement is largely free from nuclear background (the 
same assumption was made in a later reevaluation of the 
same dataset 12311 ). While the astrophysical S- factor de- 
rived in Ref. [2l| seems to agree well with theoretical 
predictions at and above the resonance, a puzzling result 
emerged below the resonance: the experimental data sug- 
gest a rather constant S-factor in the astrophysically in- 
teresting region below 400 keV; most theoretical curves, 
however, drop with decreasing energy [zjj . As we will 
show in the present paper, it is likely that this constant 
S'-factor is due to nuclear processes that cannot be ig- 
nored. 

We report in this article on a new break-up mea- 
surement performed at the SIS- 18 heavy-ion synchrotron 
at GSI (Helmholtzzentrum fur Schwerionenforschung in 
Darmstadt, Germany) using a higher-energy (150 A 
MeV) 6 Li beam. This higher beam energy should have 
several advantages compared to the one used by Kiener 
et at: (i) the stronger forward focusing allows for a more 
complete angular coverage; (ii) CD should be enhanced 
relative to the nuclear contribution. In addition, we have 
developed a comprehensive theoretical model of electro- 
magnetic and nuclear break-up processes that allows to 



interprete the measured data in detail. We will show be- 
low, however, that it is unfortunately not possible to sep- 
arate experimentally electromagnetic and nuclear contri- 
butions. Nevertheless, most of the features of the mea- 
sured data can be well explained by our model, thus giv- 
ing our calculated 2 H(a,7) 6 Li cross sections a firm ex- 
perimental basis. 



II. THEORETICAL PREDICTIONS 
A. Radiative-capture reaction 

The cross section of the 2 H(a, 7) 6 Li reaction at ener- 
gies E cm < 1 MeV is dominated by radiative E2 capture 
from d waves in the a+ 2 H channel into the J 71 " = 1 + 
ground state of 6 Li via a prominent 3 + resonance at 
E cm = 0.711 MeV. In comparison, El transitions from p 
waves to the 6 Li ground state are strongly suppressed by 
the isospin selection rule for N = Z nuclei due to the al- 
most equal charge-to-mass ratio of the deuteron and the 
a particle. Only at very low energies (E cm < 150 keV), 
the El contribution is expected to become larger than 
the E2 capture since the penetrabilities in p and d waves 
exhibit a different energy dependence [17| . 

In the past, a number of different theoretical ap- 
proaches were considered for the calculation of the low- 
energy 2 H(a,7) 6 Li capture cross section, see [H| and re- 
ferences therein. They comprise, e.g., simple potential 
models and microscopic cluster models using the resona- 
ting group method (RGM) or the generator coordinate 
method (GCM). Provided that the parameters of these 
models are well fitted to observable quantities like the 
binding and resonance energies in the 6 Li system and 
that the asymptotic form of the bound state wave func- 
tion is correctly taken into account, good agreement be- 
tween the predictions for the low-energy cross section is 
found. This shows that the radiative capture at energies 
below the 3 + resonance is essentially an extranuclear pro- 
cess and that details of the interior wave function are less 
important. 

In the present work we employ a potential model 
for the 6 Li system that provides the S'-factor for the 
2 H(a, 7) 6 Li reaction and, at the same time, can be used 
in modelling the breakup reaction when 6 Li is scattered 
on a Pb nucleus at 150 A MeV. This model assumes that 
6 Li is described by two interacting clusters, a and 2 H, 
without internal structure. Bound and scattering wave 
functions in the relevant partial waves, characterized by 
orbital angular momentum I and total angular momen- 
tum J, are obtained by solving the appropriate radial 
Schrodinger equation with a- 2 H potentials that contain 
a central part of Woods-Saxon form 
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FIG. 1: (Color online) Phase-shift data measured for low- 
energy a- 2 H scattering as a function of the relative a- 2 H en- 
ergy in the cm. system, E re i. Data points are from Jenny et 
al. (circles, [I^]), Mclntyre et al. (squares, [I(|), and Gruebler 
et al. (diamonds, 27]). The results of the model calculations 
(full lines) were obtained with the potential parameters de- 
scribed in the text. 



and a spin-orbit part of derivative Woods-Saxon form 

-l 



L-S d 



h 2 r 



dr 



1 + exp 



r-R 



(2) 



with A = 2 fm. The radius R is given by R = 
1.25 A 1 / 3 fm with A = 6; a = 0.65 fm denotes the dif- 
fuseness parameter. The depths V l c ' J were set to 60.712 
MeV for the ground state (I = 0) and 56.7 MeV for all 
other partial waves (/ = 1,2). For the spin-orbit part 
V s 'o a depth of 2.4 MeV was used for the relevant partial 
waves I = 1,2. These values were obtained by adjusting 
the parameters as to reproduce the experimental values 
for the binding energy Eb = 1.474 MeV of 6 Li and the 
3 + resonance energy with respect to the a+ 2 H threshold. 
This choice of parameters also describes the low-energy 
a- 2 H experimental scattering phase shifts very well, see 
Fig. [TJ The cross section of the radiative capture reac- 
tion is calculated in the present model with the usual 
long-wavelength approximation of the El and E2 mul- 
tipole operators M(EX^i) — er^ D Y\^(f a d) where e 
denotes the electron charge and r ad is the radius vector 
between a and deuteron. The effective charge numbers 
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depend on the charge numbers Zi and masses mi of the 
two clusters. The El effective charge number does not 
vanish since experimental values for the masses are used. 
However, the E2 contribution dominates over most of the 
range of energies with a pronounced peak at the position 
of the 3 + resonance. Only at energies below 110 keV 
the El contribution exceeds the E2 contribution. We 



display the energy dependences of the two relevant mul- 
tipole contributions to the 5-factor below in Fig. |H] of 
Sect.IV. 



B. Breakup reaction 

The theoretical description of the breakup reaction 
208 Pb( 6 Li,a 2 H) 208 Pb is considerably more involved than 
that of the radiative capture reaction, in particular if 
both electromagnetic and nuclear breakup have to be 
included. The differential breakup cross section in the 
6 Li- 208 Pb cm. system can be written in the general form 
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with reduced masses /iy = minij/frni + mj) and rela- 
tive momenta ftj = fj+j {pi/rrii — Pj/mj). fiy denotes 
the solid angle for the scattering of particles i and j in 
their cm. system and E ad = p^J {2\i ad ) is the cm. en- 
ergy in the fragment system after the breakup. In the 
initial state, 6 Li is in the ground state with total angular 
momentum Jn — 1 and Mn = ±1, 0. In the final state, 
the deuteron carries spin 1 with projections M d = ±1,0. 
The cross section (|4]) determines the relative probability 
to find the two fragments with given momenta in the fi- 
nal state und thus can be used directly in a Monte- Carlo 
simulation of the breakup reaction. 

The main task is to calculate the T-matrix element 
that contains all the relevant information on the breakup 
process. In distorted wave Born approximation (DWBA) 
it is given by 
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with the 6 Li ground state wave function <&n and the wave 
function J for the relative motion of the fragments in 
the continuum. These two functions are given by the 
solutions of the Schrodinger equation as in the calcula- 
tion of the radiative-capture cross section. The distorted 
waves x describe the scattering of the projectile on 
the target. They can be found by solving the Schrodinger 
equation for the Li-Pb scattering with the optical poten- 
tial U LiPb that only depends on the distance between Li 
and Pb. In contrast, Vi^iPb is the full many-body inter- 
action potential. It is approximated by 



Vl 



Z a Zpbe 2 



iPb 



rpb 



+U N 



Pb\ 



Z d Zpi,e 2 
\f d ~ r Pb 
rpb\) 



(6) 



U^ Pb {\r d -r Pb \) 



separating Coulomb and nuclear contributions and intro- 
ducing nuclear optical potentials U^ Pb and U^ Pb for the 
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a-Pb and 2 H-Pb interaction, respectively. Similarly we 
have 
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with rLiPb = fn — f Pb . Since both potentials contain 
Coulomb and nuclear contributions additively, it is pos- 
sible to separate the T-matrix element into a Coulomb 
and a nuclear part as 
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In general, Coulomb and nuclear contributions to the 
breakup amplitude can interfere. 

In the breakup experiment, the projectile velocity 
VLiPb relative to the target is large and the fragments 
are observed at small forward scattering angles with res- 
pect to the beam axis. Thus it is sufficient to replace 
the distorted waves appearing in Eq. §5§ by their eikonal 
approximation, i.e. 
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where the coordinate vector has been decomposed as 

rLiPb — b + ze b0 am b _L e bca m (12) 

and ebeam denotes the beam direction. The Coulomb 
part of the phase function can be calculated analytically. 
In order to avoid a divergent result at small impact pa- 
rameters, b, the Coulomb potential of a point-like target 
charge in Eq. ([7]) is replaced by that of a homogeneous 
sphere. In the eikonal approximation, the Coulomb and 
nuclear T-matrix elements can be written as 
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with the Coulomb form factor 
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that can be both decomposed into multipoles L = 

0, 1,2, Using the method of steepest descent, the mul- 

tipole components of the Coulomb form factor (fT4|) are 
easily calculated. Neglecting the nuclear contribution in 
the phase function SLiPb one obtains in lowest order the 
well-known Coulomb excitation functions in the semiclas- 
sical approximation. The nuclear optical potential U^ iPb 
in (ITT1) leads to corrections that take, e.g., the absorp- 
tion by the target nucleus into account; also, relativistic 
corrections are easily included, see Ref. [28| for details. 

Since the E2 virtual photons are orders of magnitude 
more abundant than the El ones, the Coulomb contri- 
bution to the breakup is essentially sensitive only to the 
quadrupole contribution. (There is no monopole con- 
tribution in this case.) Contrary to electromagnetic El 
excitations, nuclear L — 1 excitations are not suppressed 
by isospin selection rules. Nuclear processes for all multi- 
polarities have, therefore, to be taken into account when 
modelling the break-up of 6 Li into a+ 2 H at about 150 A 
MeV. In the present work, we included nuclear L — 0, 1, 2 
excitations because higher multipoles are expected to 
give only small contributions to the total breakup am- 
plitude. 

In order to obtain numerical results for the T-matrix 
elements, one has to specify the nuclear optical potentials 
that enter into the calculation. Unfortunately, there are 
no systematic optical- model potentials available descri- 
bing the elastic scattering of a, 2 H and 6 Li on a Pb tar- 
get at 150 A MeV. Therefore, we generated the optical 
potentials from systematic optical-model potentials for 
nucleon-Pb elastic scattering and folded them with the 
matter distribution of the projectile and the fragments, 
respectively. These potentials were tuned to reproduce 
published elastic-scattering data at incident energies as 
close to 150 A MeV as possible by multiplying the real 
and imaginary parts by scaling factors not too far from 
unity. Literature data have been used for the elastic scat- 
tering of 2 H + 208 Pb at 55 and 70 A MeVH HI. of a 
+ 208 Pb at 120 A MeV and 175 A MeV [ijand of 6 Li 
+ 208 Pb at 100 A MeV (H. 



We found that deuteron and a scattering on Pb were 
best described starting with the relativistic nucleon- 
nucleus potentials of Ref. 33]. In the case of 6 Li- 
Pb scattering the non-relativistic optical-model poten- 
tial from Ref. [lij for nucleon-nucleus scattering worked 
best. Fig. [5] shows measured and fitted elastic-scattering 
data for the three cases. The optical-model potentials, 
obtained by the procedure described above for the actual 
energy of the breakup experiment, are well fitted by a 
Woods-Saxon shape. Since mostly the outer region of 
the potential is important the fits were started at a ra- 
dius of 7 fm. In Table U we give the numerical values of 
the depth, radius and diffuseness parameters for the real 
and imaginary parts. 
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TABLE I: Woods-Saxon potential parameter used to describe 
the scattering of 6 Li, a, and 2 H on a Pb target. 
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FIG. 3: (Color online) Expected distribution of the differen- 
tial cross section, da /cL9q, as a function of the scattering angle, 
#6, of the excited 6 Li* before break-up, in arbitrary units. The 
full (red) curve represents the total distribution, whereas the 
nuclear and Coulomb contributions are depicted by the dot- 
dashed (blue) and dashed (green) histograms, respectively. 
Note that the different curves have been normalized to the 
same total cross section. All distributions were summed over 
2 H-a cm. energies, E re i, up to 1.5 MeV. 



FIG. 2: (Color online) Center-of-mass angular distributions 
for (a) 55 and 70 A MeV 2 H, (b) 120 and 175 A MeV a, 
and (c) 100 A MeV 6 Li on 208 Pb. The full lines represent fits 
to the measured data using the optical-model potentials as 
described in the text. Note that the angular distributions for 
E d = 140 MeV in panel (a) and for E a = 699 MeV in panel 
(b) have been scaled by a factor of 10~ 2 . 



C. Predicted observables 

The most meaningful observable that can illustrate the 
predictions from the above-sketched model of 6 Li break- 
up is the scattering angle, 9q, of the excited 6 Li* before 
break-up, relative to the incoming 6 Li beam. Fig. [3] de- 
picts the expected 8$ distribution. 

The figure clearly shows that pure nuclear, pure 



Coulomb, and total (CD+nuclear) distributions exhibit 
distinctly different peak structures. Pure Coulomb in- 
teraction has its most prominent peak where the other 
contributions show a minimum. Likewise, the total 
(CD+nuclear) distribution can be distinguished from a 
nuclear-only theory by the large amplitude of the most 
prominent peak (due to constructive CD-nuclear interfe- 
rence), and by the disappearance of the third maximum 
(due to destructive interference). In principle, these fea- 
tures should allow to separate the contributions from the 
individual interactions. However, the theoretical predic- 
tions have to be folded with the resolution and the ac- 
ceptance of the experimental apparatus using the Monte- 
Carlo simulations described below in subsection IIIIB1 
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III. EXPERIMENTAL PROCEDURE 

A. Apparatus 

A schematic view of the set-up used is given in Fig. |4] 
A 208 Pb target with 200 mg/cm 2 thickness was bom- 
barded by a primary 6 Li beam of 150 A MeV energy. 
The 6 Li beam was produced by the SIS-18 synchrotron 
at GSI, separated from possible contaminant ions by us- 
ing the FRS FRagment Separator [35[ and transported 
to the standard target position of the kaon spectrometer 
KaoS [3fj. The average 6 Li beam intensity at the break- 
up target was of the order of 5xl0 4 per 4 sec spill. The 
beam had a width of 0.17(0.12) cm and an angular di- 
vergence of 4.4(4.4) mrad in x(y) direction at the target 
(la widths). 

The angles and positions as well as the energy losses of 
the outgoing particles, 2 H and a, were measured by two 
pairs of single-sided silicon strip detectors (SSD, 300 
thick, 100 /im pitch) [33] placed at distances of 15 and 
30 cm, respectively, downstream from the target. From 
the detector pitch one can calculate a resolution of the 
2 H-a opening angle in the laboratory, #24, of about 1%. 
Non-interacting 6 Li beam particles were identified event 
by event with a 16-strip AE detector located directly be- 
hind the SSD and stopped in a cylindrical Ta absorber 
(12 mm diameter, 20 mm length) placed behind the de- 
tector. Break-up events were discriminated from non- 
interacting 6 Li-beam events by their energy-loss signals 
in the 16-strip AE detector; an energy loss correspond- 
ing to 6 Li was used as a trigger veto signal. Deuteron 
and a momenta were analyzed with the large-acceptance 
KaoS spectrometer and were detected in two consecutive 
multi-wire proportional chambers (MWPC [37]]) followed 
by a plastic-scintillator TOF wall consisting of 30 ele- 
ments (each 7 cm wide and 2 cm thick). This plastic 
wall was used as a trigger detector for the data acquisi- 
tion system. The KaoS magnets' volume was filled with 
He gas at atmospheric pressure to reduce multiple scat- 
tering. 

The coincident 2 H and a fragments resulting from 
break-up in the 208 Pb target were identified by recon- 
structing their vertex at the target. This removed all 
break-up events in layers of matter other than the target. 
The 2 H and a momenta were determined from tracking 
them with GEANT through the MWPC and TOF wall 
behind KaoS. The incident angles in front of the magnets 
were known unambiguously from the SSD hits. While in 
the SSD each hit could be attributed to either 2 H or 
a by its energy deposition, the corresponding hits in the 
MWPC were attributed to the respective particle type by 
finding the optimum trajectory through the MWPC and 
the TOF wall. This was done in an iterative procedure 
that started with a test assignment of each hit to cither 
a or 2 H and a test momentum for each of them. Both 
the momentum values and the assignments were then it- 
eratively changed until the minimum squared deviation 
from the observed hits in all detectors downstream from 



the KaoS magnet were reached. This momentum recon- 
struction could be shown to be accurate within about 
10~ 3 . From the opening angles between the fragments 
and from their momenta, the relative energies, E re i, be- 
tween the 2 H and a particles in the cm. system could 
be reconstructed. 



B. Monte-Carlo simulations 

It is obvious that the experimental apparatus imposes 
strong restrictions on the detection of the break-up par- 
ticles, a and 2 H. This applies in particular to the angular 
acceptance, the energy and position resolution, and the 
detection efficiency. As a consequence, a meaningful com- 
parison between theoretical predictions and experimental 
data can only be made using theoretical data filtered by 
the experimental set-up. To this end, we have modelled 
the entire set-up, starting in front of the 208 Pb break-up 
target, in GEANT3 38]. As an event generator, the the- 
oretical model described in the previous section was used. 
Input data were generated as statistically-distributed en- 
sembles of 100,000 break-up "events" each that were dis- 
tributed according to the calculated differential cross sec- 
tions. The emittance of the 6 Li beam (as measured with- 
out break-up target and without absorber) was imposed. 
Each break-up particle, a and 2 H, was followed through 
the remainder of the Pb target after the reaction vertex, 
the SSD detectors, the beamline exit window, the He- 
filled interior of the magnets with the magnetic field and 
the air behind KaoS before hitting the MWPC volumes 
and the TOF wall. 

The Monte-Carlo simulations were used to obtain esti- 
mates of the resolution and the efficiency of our setup. As 
an example, we plot in FigfSJa) the ler-resolution of E re i. 
The data points were obtained by sending 10,000 events 
each with different values of E re i (within a narrow bin 
of 0.1 MeV width for each case) into our setup and an- 
alyzing the outgoing particles with the same routines as 
in the experiment. From the same data sets, the number 
of counts gave an approximate estimate of the detection 
efficiency, shown in Fig(5jb). In the experiment, how- 
ever, the detection efficiency is additionally limited by 
the small and strongly fluctuating energy deposition of 
deuterons in the MWPC. This latter quantity cannot be 
simulated easily, such that we had to normalize the num- 
ber of observed and simulated counts. Therefore, our 
experiment does not allow to determine absolute cross 
sections, despite the fact that all incident 6 Li ions were 
counted. 



IV. EXPERIMENTAL RESULTS 

A relatively unbiased observable, based only on high- 
resolution SSD measurements, is the opening angle, 6*24, 
between the outgoing fragments 2 H and a. Fig. [5] shows 
this distribution, summed over E re i values up to 1.5 
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FIG. 4: (Color online) The experimental setup shows the fragment-tracking SSD behind the Coulomb-break-up target followed 
by a 16-strip AE detector and a beam stopper. Deuteron and a positions were measured near the focal plane of the KaoS 
QD-spectrometer by two successive large area multi-wire proportional chambers (MWPC) followed by a scintillator-paddle 
TOF wall used for trigger purposes. 
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FIG. 5: (a) Relative-energy resolution (la widths) as deter- 
mined by simulating 6 Li breakup with GEANT. The data 
points can be approximated by the fitted function oe = 
0.1645 x y/E re i. (b) Combined geometrical and analysis ef- 
ficiency of determining E re i from the 2 H and a momentum 
vectors. The intrinsic efficiency of the MWPC detectors has 
been assumed to be unity in this graph. 



MeV; this condition was also set for all other spectra 
shown below. The experimental data points are com- 
pared with the corresponding Monte-Carlo simulations 
for pure Coulomb (CD) and pure nuclear interactions as 
well as combined (CD+nuclear) interactions. Each sim- 
ulated histogram was normalized to contain the same 
number of counts as the experimental spectrum, thus 
providing the single scaling factor used to normalize all 



simulated distributions. 
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FIG. 6: (Color online) Opening angles, 824, between the out- 
going fragments 2 H and a. Full circles correspond to mea- 
sured data. The dash-dotted histogram (blue) denotes sim- 
ulations with pure nuclear interaction, whereas the pure CD 
contribution is shown by the dashed histogram (green). Com- 
bined (CD+nuclear) contributions are shown by the full red 
line. Note that the numbers of simulated counts in each spec- 
trum were normalized to the experimental ones. 

As seen in Fig. [BJ the data are reasonably well repro- 
duced by the simulations over their entire range. The 
3 + resonance peak is clearly visible around 3 degrees; its 
angular width is well reproduced indicating that the sim- 
ulation takes both the scattering and the finite angular 
resolution well into account. It is obvious, however, that 
one cannot distinguish between the different interactions 
on the basis of this angular distribution. We have there- 
fore to search for an observable that is more sensitive to 
the type of interaction. In Fig. [3] above, we have shown 
that the observable 9$ should be very sensitive to the 



8 



40 

30 

20 

10 


150 



3 100 




(a) 

• Experiment 
— CD+Nuclear - 

■■ CD 
■ — Nuclear 



(b) - 



r 1 *l 



r j I i i jr iTH 




FIG. 7: (Color online) Angular distribution of the excited 
6 Li* nuclei after the reaction {9$). The panels represent 
three different bins of E rel : (a) 0.0 < E rel < 0.5 MeV; (b) 
0.5 < E re i < 0.9 MeV; (c) 0.9 < E rei < 1.5 MeV. The mea- 
sured data points are shown in comparison with simulations 
with pure nuclear and pure CD as well as with (CD+nuclear) 
interactions. Line types and color codes are identical to the 
ones in Fig. [6] 



type of interaction. 

The experimental data for this observable are pre- 
sented in Fig. [7J Panel (a) of Fig. [JJ shows data for 
E re i values below the resonance; panel (b) covers the 
resonance region, whereas panel (c) has been plotted for 
0.9 < E re i < 1.5 MeV. The finite KaoS acceptance cuts 
the distributions at about 4 degrees. The figure shows 
clearly that the observable 9$ is sensitive to the type of 



interaction. In all panels, the combined (CD+nuclear) in- 
teraction, including interference, reproduces most of the 
structures observed in the data points (red histograms). 
This is particularly true for the sub-resonance region, 
panel (a). The green histograms (CD-only) show sin- 
gle peaks at larger angles. The pure nuclear interaction 
(blue histograms) rises rapidly at small values of 9q, in 
agreement with the measured data, but lacks the struc- 
tures visible in the data points. The narrow peaks visible 
in the on-resonance data, panel (b), at values of 06 of 
fa 1.5,2.6 and 3.3 degrees, are not perfectly reproduced 
by the (CD+nuclear) model and point to small deficien- 
cies of the theoretical model. Nevertheless, Fig.UJdemon- 
strates that Coulomb-nuclear interference is at work and 
that the signs of the interference terms are correct. We 
also conclude that even at our incident energy of 150 A 
MeV the nuclear break-up is dominant. 

The angle-integrated energy-differential cross sections 

as a function nf K„„7 a.rp shown in Fie |^] Thp full his- 



800 




E re| (MeV) 

FIG. 8: (Color online) Differential cross sections as a function 
of the energy, E re i, in the a- 2 H cm. system. Black points 
indicate the experimental data; the histogram corresponds to 
the GEANT simulation using the (CD+nuclear) interaction 
as described in the text and a binning of 100 keV (note that 
the vertical error bars result from a quadratic sum of statis- 
tical and systematical uncertainties). 

togram was obtained from the (CD+nuclear) calculation 
convoluted by our GEANT simulation and normalised to 
the experimental yield. The points and the histogram 
represent the measured and predicted differential cross 
sections, respectively, as a function of E Te \. Our E re i 
distribution is in very good agreement with the simula- 
tion in particular in the energy region below 400 keV. 
As we will show below (see Fig fTOj) , the differential cross 
sections in this energy regime result mostly from nuclear 
interactions. 

The astrophysically important quantity is the astro- 
physical iS-factor, 524, for the 2 H(a, 7) 6 Lireaction. Since 
nuclear processes dominate, in particular for low E re i, 
the determination of this quantity via an evaluation of 
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FIG. 9: (Color online) (a) Theoretical El, E2, and total S 2 4- 
factors that describe well the present experimental data, to- 
gether with data points from the previous CD experiment by 
Kiener et al. [2lJ and from direct measurements (Robertson 
et al. fjjj and Mohr et al. 03). See Sect |Vj for an inter- 
pretation of the data of Ref. [21j. (b) Comparison of various 
theoretical predictions for the summed El- and E2 contribu- 
tions to S 2 4(E) QjHHIMS 



the CD component in our break-up data is not feasible 
since. However, we have demonstrated above that our 
theoretical model describes well the measured cross sec- 
tions, hence the astrophysical S-factors from the present 
work are those from our theoretical model. 

The resulting E1-, E2-, and total 5'24-factors are vi- 
sualized in the upper part of Fi g. E] together with the 
previous CD data of Kiener et al. [21| and the direct data 
of Mohr et al. [IH and Robertson et al. [13] ■ The present 
results for the £ , 2-component are in good agreement with 
the direct measurements of Refs. fijl Ral in the resonance 



region and above which gives confidence in our model. 

Another check of the validity of our treatment of the 
Coulomb part of the break-up reaction, described in 
Sect lll Al can be done by comparing our calculated res- 
onance parameters for the 3 + resonance with the exper- 
imental ones. In order to determine those parameters, 
we have calculated the theoretical capture cross section 
around the resonance in 1 keV steps and fitted a Brcit- 
Wigner parametrization to the resonance. We obtain T- 
widths of T a = 22.1 keV and T 7 = 0.437 meV, in good 
agreement with the literature values of T a = 24 ± 2 keV 
and T 7 = 0.440 ± 0.030 meV as cited by Mohr et al. [H]. 
Note that we have used a spectroscopic factor of unity. 
We will comment on the data points from the previous 
CD experiment [2l[ in the following section. 

The direct 2 H(a, 7) 6 Lireaction at very low energies is 
sensitive also to the El amplitude. In our experiment, 
this component cannot be constrained experimentally 
due to the weak flux of virtual El photons. We have 
therefore to rely exclusively on the theoretical model. At 
higher energies, however, Robertson et al. [17] could sep- 
arate El- and _E2-components on the basis of measured 
angular distributions. Fig.l in [TtJ shows that their the- 
ory seems to overestimate the El component. Our El 
curve is very close to Robertson's so that we also seem 
to overestimate this component. 

Several theoretical models for 6 Li have been pro- 
posed to determine the shape and the magnitude of 
the S24 energy dependence, such as potential models 
[IH |H, [3^ , cluster- model calculations [24{ , or ab-initio 
calculations Those predicted curves for S24 which 

include both El- and _B2-contributions are displayed to- 
gether with the theoretical curve from this work in the 
lower panel of Fig. |9] As one can see in this figure, all 
the calculations shown -independent of their very diffe- 
rent model assumptions- yield very similar curves. We 
have not included the theory of Blokhintsev et al. (42J) 
because it was specifically tuned to approach the exper- 
imental data of Ref. [2l| ■ 



V. COMPARISON WITH OTHER CD 
EXPERIMENTS 

As visible in the upper panel of Fig. |9l the low-energy 
data points derived for the S^-factor from the work of 
Kiener et al. [2l[ disagree with the theoretical curve that 
we have deduced from the present work. We believe that 
this is most likely due to a strong nuclear contribution 
at the lower incident energy of 26 A MeV, which was not 
considered in their analysis in view of the theoretical pre- 
diction [2^. We have performed a calculation with the 
theoretical model of this work at an incident energy of 
26 A MeV and have calculated CD and nuclear cross sec- 
tions for the laboratory- angular range between 1.5 and 6 
degrees, which should correspond approximately to the 
acceptance of the setup of Ref. [2l|. Fig. fTTJI displays the 
ratio of nuclear to Coulomb break-up cross sections as a 
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function of E re i for the two bombarding energies. Our 
theory predicts that at 150 A MeV the nuclear cross sec- 
tions are about a factor of 3 larger than the CD ones at 
and above the resonance, whereas the nuclear component 
dominates strongly at the lowest energies. At 26 A MeV, 
the ratio of nuclear to CD cross sections is predicted to 
be about a factor of ten larger than at 150 A MeV over 
the entire range of E re i. This suggests that the data 
points shown in Ref. [2l[ result almost exclusively from 
nuclear interactions, contrary to the assumptions under- 
lying their analysis. It is therefore not very meaningful to 
tune theoretical models in order to improve their agree- 
ment with the 26 A MeV data as was done in Ref. 1431 . 
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FIG. 10: (Color online) Ratio of nuclear and Coulomb differ- 
ential dissociation cross sections for 6 Li at 150 A MeV (full 
line) and at 26 A MeV (dashed line). Both curves were calcu- 
lated with the same model described in detail above in section 
II.B. 



VI. PRODUCTION OF 6 LI IN THE BIG BANG 

The 2 H(cv,7) 6 Li reaction is the main path for 6 Li 
BBN production while destruction proceeds via the 
6 Li(p,a) 3 He one. Both rates are available in the 
NACRE [43J compilation. While the latter reaction rate 
is reasonably well known at BBN energies, prior to this 
experiment the former suffered from a large uncertainty. 
This was mainly due to the fact that the published error 
margins were aimed at including the Kiener et al. [2l[ 
measurement [44j. As a result from the present study, 
we can now propose a more reliable central value based 
on a successful theoretical model, and a safe upper limit 
that is even somewhat smaller than the previous NACRE 
upper limit. 

In Fig. [Til the BBN abundances of 7 Li and 6 Li are dis- 
played as a function of the baryonic density. (It is usual 
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FIG. 11: (Color online) Predicted BBN production ratios for 
6,7 Li over hydrogen as a function of 77, the baryon-to-photon 
ratio in the early Universe. The solid red line represents the 
result for 6 Li from the S'24-values obtained in the present 
work, based on theoretical values for the El and E2 com- 
ponents. The dashed red line represents a very conservative 
but safe upper limit where all observed events are assumed 
to result from Coulomb break-up. The blue band denotes 
the range of predicted 7 Li yields Q. Observational data are 
indicated by horizontal green-hatched areas: the upper one 
has been derived from the recent review of lithium observa- 
tions by Spite and Spite [f| ; the lower one corresponds to the 
largest 6 Li yield reported for the star HD 84937 O. The 
yellow vertical band shows the WMAP 77- value [lj]. 



to introduce another parameter, rj, the ratio of the num- 
ber of baryons over the number of photons which remains 
constant during the expansion, and which is directly re- 
lated to fl b by f7 fc /i 2 =3.65xl0 7 ?7). The blue 7 Li band is 
the result of a Monte Carlo calculation taking into ac- 
count nuclear uncertainties as described in Ref. Q. The 
upper hatched horizontal area in the figure shows the pri- 
mordial lithium ( 6 Li+ 7 Li) abundance derived from the 
"Spite plateau", i.e. from the practically metallicity- 
independent Li observations in metal-poor stars 

The solid red line for 6 Li has been calculated within the 
same physical model, using the theoretical reaction rate 
from this work for 2 H(a,7) 6 Li. Both El and E2 contri- 
butions have been included. At WMAP baryonic density, 
a value for the 6 Li/H production ratio of s» 1.3 x 10~ 14 
results. The dashed red line represents a very conserva- 
tive upper limit for 2 H(a, 7) 6 Li that would hold if the 
low-energy S24 data points from this work would result 
from CD only. Fig. [7] demonstrates clearly that this is 
not the case. But even this extremely conservative limit 
is about two orders of magnitude smaller than the only 
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TABLE II: Recommended 2 H(a, 7) 6 Li reaction rates using 
theoretical El— and E2 — S-factors from the present work for 
the temperature range 10 6 K < T < 10 10 K (10 -3 < T 9 < 10). 

positive observation of 6 Li surviving after the reanalysis 
of Li lines by Steffen et al. [X3j | . indicated in Fig. ITT1 by 



the lower hatched horizontal band. This finding corrobo- 
rates earlier statements (e.g. [3, that observations -if 
confirmed- of 6 Li primordial yields around a few percent 
of the Spite plateau would require astrophysical sources 
other than BBN. 

In order to facilitate astrophysical calculations of stel- 
lar 6 Li synthesis with our new theoretical El and E2 
S- factors, we list in Table Ql] the reaction rates for the 
temperature range 10 6 K <T< 10 10 K. 

VII. CONCLUSIONS 

A kinematically complete measurement of the high- 
energy break-up of 6 Li at 150 A MeV has shown that 
Coulomb and nuclear contributions and their interferen- 
ces have to be taken into account when interpreting the 
measured angular distributions. Though it was not possi- 
ble to extract the Coulomb part experimentally, we were 
able to infer the E2 component of the astrophysical S24- 
factor for the 2 H(a,7) 6 Li reaction from a theoretical re- 
action model that describes well in particular the low- 
energy break-up data. The model predicts a drop of S24 
with decreasing relative 2 H-a energy, E re i, as predicted 
also by most other nuclear models for 6 Li, contrary to 
conclusions from an earlier CD experiment performed at 
the lower energy of 26 A MeV. We have presented evi- 
dence that this earlier experiment probably has measured 
mostly nuclear break-up of 6 Li. Our findings allow to 
make new predictions for the 6 Li/H production ratio in 
Big-Bang nucleosynthesis (BBN) which is orders of mag- 
nitudes smaller than the one derived from claimed obser- 
vations of 6 Li in old metal-poor stars. Sources other than 
BBN have therefore to be invoked for 6 Li production if 
those observations are confirmed. 
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